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The functional RG for the random field and random anisotropy 0{N) sigma-models is studied to two loop. 
The ferromagnetic/disordered (F/D) transition fixed point is found to next order in d = 4 + e for A'^ > A'^^ 
{Nc = 2.8347408 for random field, Nc = 9.44121 for random anisotropy). For N < Nc the lower critical 
dimension d = die plunges below die = 4: we find two fixed points, one describing the quasi-ordered phase, 
the other is novel and describes the F/D transition, die can be obtained in an {Nc — A'^)-expansion. The theory 
is also analyzed at large and a glassy regime is found. 



It is important for numerous experiments to understand 
how the spontaneous ordering in a pure system is changed by 
quenched substrate impurities. One class of systems are mod- 
eled by elastic objects in random potentials (so-called ran- 
dom manifolds, RM). Another class are 0{N) classical spin 
models with ferromagnetic (Ferro) couplings in presence of 
random fields (RP) or anisotropies (RA). The latter describe 
amorphous magnets QJ. Examples of RF are liquid crystals 
in porous media f^, He-3 in aerogels |3|, nematic elastomers 
(a], and ferroelectrics |5]. The XY random field case N = 2 
is common to both classes and describes periodic RM such as 
charge density waves, Wigner crystals and vortex lattices |6]. 
Larkin showed |7] that the well-understood pure fixed points 
(FP) of both classes are perturbatively unstable to weak disor- 
der for d < due (due = 4 in the generic case). For a contin- 
uous symmetry (i.e. the RF Heisenberg model) it was proven 
10] that order is destroyed below d = 4. This does not settle 
the difficult question of the lower critical dimension d\c as a 
weak-disorder phase can survive below due, if associated to 
a non-trivial FP, as predicted in d = 3 forjhe Bragg-glass 
phase with quasi long-range order (QLRO) |9]. For the ran- 
dom field Ising model = 1 (RFIM) it was argued 1 10|, then 
proven fllll that the ferromagnetic phase survives in d — 3. 
Developing a field theory to predict die, and the exponents of 
the weak-disorder phase and the Ferro/Disordered (F/D) tran- 
sition, has been a long-standing challenge. Both extensive nu- 
merics and experiments have not yet produced an unambigu- 
ous picture. Among the debated issues are the critical region 
of the 3D RFIM JlJ and thepossibility of a QLRO phase in 



A peculiar property shared by both classes is that observ- 
ables are identical to all orders to the corresponding ones in a 
d — 2 thermal model |^5|. This dimensional reduction (DR) 
naively predicts die = 4 for the weak-disorder phase in a 
RF with a continuum symmetry fld\ and no Ferro order for 
the d = 3 RFIM, which is proven wrong |11]. It also pre- 
dicts d^e = 6 for the F/D transition FP. While there is agree- 
ment that multiple local minima are responsible for DR fail- 
ure, constructing the field theory beyond DR is a formidable 
challenge. Recent attempts include a reexamination of the (j)^ 
theory (i.e. soft spins) for the F/D transition near d = 6 fl7i. 



Previous large- approaches failed to find a non-trivial FP, 
but a self-consistent resummation including the 1/N correc- 
tions hinted at exponents different from DR (without succeed- 
ing in computing them) from a solution breaking replica sym- 
metry 1 18]. 

As for the pure 0{N) model, an alternative to the soft-spin 
version (near d = 6) is the sigma model near the lower critical 
dimension (here presumed to be d = 4). In 1985 D.S. Fisher 
lligll noticed that an infinite set of operators become relevant 
near d = 4 in the RF 0{N) model. These were encoded 
in a single function for which Functional RG equations 
(FRG) were derived to one loop, but no new FP was found. 
For a RM problem fid] it was found that a cusp develops in 
the function (the disorder correlator), a crucial feature 
which allows to obtain non-trivial exponents and evade DR. 
A fixed point for the RF model was later found ]9] in d = 
4 — e for = 2. It was noticed only very recently ]21] that 
the 1-loop FRG equations of Ref. 111911 possess fixed points in 
d = 4 + e for iV > 3, providing a description of the long- 
sought critical exponents of the F/D transition. 

In spite of these advances, many questions remain. Con- 
structing FRG beyond one loop (and checking its internal con- 
sistency) is highly non-trivial. Progress was made for RM 
1I22I l23ll . and one hopes for extension to RF. Some questions 
necessitate a 2-loop treatment, e.g. for the depinning transi- 
tion, as shown in L24J . In RF and RA models the 1 -loop anal- 
ysis predicted some repulsive FP in d = 4 + e (for larger val- 
ues of N), and some attractive ones f9','2i| in d = 4 — e. The 
overall picture thus suggests a lowering of the critical dimen- 
sion, but how it occurs remains unclear Finally the situation 
at large N is also puzzling. Recently, via a truncation of exact 
RG ]26] it was claimed that DR is recovered for N large. 

Our aim in this Letter is twofold. We reexamine the overall 
scenario for the fixed points and phases of the 0{N) model 
using FRG. This requires the FRG to two loop. Here we 
present selected results, details are presented elsewhere ]27^. 
We find a novel mechanism for how the lower critical dimen- 
sion is decreased below d = 4 for N < Nc at some critical 
value Nc- We obtain a description of the bifurcation which 
occurs at Nc, and below Nc we find two perturbative FPs. 
Thanks to 2-loop terms die can be computed in an expansion 
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in Nc — N, and the Ferro/Para FP below d = 4 is found. A 
study at large N indicates that some glassy behavior survives 
there. 

Let us consider 0{N) classical spins n{x) of unit norm 
= 1. To describe disorder-averaged correlations one in- 
troduces replicas na{x), a = 1, . . . ,k, the limit fc = being 
implicit everywhere. The starting model is a non-linear sigma 
model, of partition function Z = f V^k] and action: 



To 



(1) 



ab 



where Ha = {<Ja,'^a) with aa{x) — \/l ~ 7fa(x)2. A small 
uniform external field ^ Mo(l, 0) acts as an infrared cutoff. 
Fluctuations around its direction are parameterized by (A^ — 1) 
TT-modes. The ferromagnetic exchange produces the 1 -replica 
part, while the random field yields the 2-replica term Ro{z) ~ 
z for a bare Gaussian RF. RA corresponds to Ro(z) — z^ . As 
shown in 1 19] one must include a full function Ro{z), as it is 
generated under RG. It is marginal in d = 4. 

To obtain physics at large scales, one computes perturba- 
tively the effective action r[na (x)] . It can be expanded in gra- 



dients near a uniform background configuration n", and split- 
ted in 1-, 2- and higher-replica terms. From rotational invari- 
ance it is natural to look for T in the form Q with Ha ^ = 

(df , Tf^), aa ^ ^ - W)^ - TTf = Z" V^^,, 

To^Tb. = Tq/Zt, Mo Mr = MqV^/^t, m = ^/M^ 
the renormalized mass of the 7?a modes, and Ro{naHb) 
nfR{Ha H^). Higher vertices generated under RG are irrele- 
vant by power-counting, hence discarded. Renormalization of 
T contributes to the flow of R, and one sets T = at the end. 

One computes Z, Zt and R perturbatively in Ro and 
extracts (3 and 7 functions /3[i?](z) = —mdmR{z), 7 = 
—nidm In Z and jt = —mdm In Zt, derivatives taken at 
fixed Ro,To, Mo- Although calculation of the Z-factors is 
simplified due to DR, anomalous contributions appear from 
the non-analyticity of R{z). To compute R{z), one chooses 
a pair of uniform background fields n°) for each (a, b). 
We use a basis for the fluctuating fields (to be integrated over) 
such that Ha = {aa,ria,Pa)' = (cfo, rib, pb), where 77 Hes in 
the plane common to (n" , n°), and pa along the perpendicular 
N — 2 directions; both have diagonal propagators. Denoting 



cos (j)ab, one has HaHb = cos (j)ab {(JaCrb + IhVb) 



sin^afc (cTaiib - (Jbrja) + Papb- One gets factors of [N - 2) 
from the contraction of p. Our calculation to 2 loops results in 
the flow-equation for the function i?(0) = R{z — cos 0), and 

e = 4 - d: 



dtR{(j>) 



eR{4>) 
1 



'-R"{(j>f - R"{0)R"{4>) + (iV-2) 



{R"{(t)) - R"{Q))R"'[(^f + (iV-2) 



cot ( 



1 R'{<i>Y 



-R'{4>? 



-cot (f>R'{(f>)R"{0) 
5 + cos 2(j) 



R'{<I>?R"W 



2sin^ 



-R"i<P)' 



4 sin* I 



N+2 



R"'{0+fR"{c^) 



■i?"(0)(2(2 + cos2(^)i?'((/))2 

N 



4 sin* , 

Qsm2(j)R'{(j))R"{(j)) + (5 + cos 20) sin^ (pR!'{<P)' 



cot (j)R"'{Q+YR'{(l)) - 2{N-2) R"{0) - R"{Of + 7„i?"'(0+)2 R{^) (2) 



with di := —mdm, and the last factor proportional to R{(f>) is 
— 27T and takes into account the renormalization of temper- 
ature. Thanks to the anomalous terms, arising from a non- 
analytic R{(t>), this /3-function preserves a (at most) linear 
cusp (i.e. finite R"'{{)^)), and reproduces for N = 2 the pre- 
vious 2-loop results for the periodic RM 123]. For > 2, 
anomalous contributions are determined following i23l . 7 is 
found as 



j^{N~ l)i?"(0) + 



3N-2 7J"Vn+^2 



8 



R"'io+y 



(3) 



either via a calculation of (aa) L3(X] . or of the mass corrections, 
a result consistent with the /? function (|2j Isill . The determi- 
nation of 7t is more delicate fs?], and we have allowed for an 
amalous contribution 7a, whose effect is minor, and discussed 
below. The correlation exponents (standard definition j^) 
are obtained as 77 = e — 7, 77 = 7t — 7 at the FP. (|3 has the 
form: 

deR^eR + 'B{R,R) + CiR,R,R) + OiR^) (4) 



r 



We now discuss its solution, first in the RF case, and setting 
7a = 0. The 1-loop flow-equation (setting C = 0) admits, 
in dimensions larger than 4, a fixed point R^^j^ with a single 
repulsive direction, argued by Feldman to describe the F/D 
zero temperature transition. This is true only for N > Nc- For 
N < Nc this fixed point disappears and instead an attractive 
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FIG. 1: (Color online) Phase diagram. D = disordered, F = ferro- 
magnetic, QLRO = quasi long-range order. 
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FIG. 2: (Color online) Parametric plot for solutions of ^5} for < 
Nc (dots) for RF, equivalent to J6j (solid line, parabola) and flow 
(arrows). / parameterizes disorder and only / > is physical. Cmp. 
with fig.Q right. 

fixed point Rqi^^q appears which describes the Bragg glass 
for = 2. We have determined = 2.8347408 and the 
solution Rc{u) which satisfies B{Rc, Rc)\N=Na = 0. It is 
formally the solution at e = 0. Since the FRG flow vanishes 
to one loop along the direction of Rc, examination of the 2- 
loop terms is needed to understand what happens at = Nc- 
In particular the F/D transition should still exist for N < Nc, 
though it cannot be found at one loop. It is not even clear a 
priori whether it remains perturbative. 

The scenario found is perturbative, accessible within a dou- 
ble expansion in y/\e\ and N — Nc- To this aim, we write the 
leading terms in — iVc and e of 0, namely 

deR = eR + B^R, R) + Cc{R, R, R) + (A^-A^c) Bjv(i?, R) 
Bc{...)^B{...)\n=n^, Cci. . .) ^ C{. . .)\n=n^ (5) 

One looks for a fixed-point solution of the form R{u) = 
gRc{u) + g^5R{u), with g > 0, i?^'(0) = -1, and its flow. 
This analysis is done numerically and leads to the flow shown 
schematically on Fig.H The RG-flow projected onto the di- 
rection of g is equivalent to 

dig = eg+ 1.092(Ar - Nc)g^ + 2.352g^ . (6) 

As a solution of the functional flow near Nc, its simplicity is 
surprising. Setting g = {Nc — N)f, there are three FP: 

^^^^-^^- 1.092/ + 2.352/2 = 0, or / = 0. (7) 

For N > Nc the physical branch is / < 0. As seen in Fig. [2 
for d > 4 there is a ferro phase (i.e. / = is attractive) and an 
unstable FP describing the F/D transition, given by the nega- 
tive branch of Q. At A^ = A'c one sees from (6) that the F/D 
fixed point is still perturbative, but in a ^/e expansion for g 
(and for the critical exponents). For N < Nc the physical side 
is / > and there are two branches on Fig.Qcorresponding 
to two non-trivial fixed points. One is the infrared attractive 
FP for weak disorder which describes the Quasi-Ordered fer- 
romagnetic phase; the second one is unstable and describes 



the transition to the disordered phase with a flow to strong 
coupling. These two fixed points exist only for e < Ec and an- 
nihilate at Cc. The lower critical dimension of the RF-model 
for N < Ncis lowered from d = 4 to 

dJJF = 4-ec w 4:-0.1268{N -Nc)^ + 0{{N-Ncf) . (8) 

Note that the mechanism is different from the more conven- 
tional criterion d — 4 + ri{d) = at d = djc. 

The same analysis for the random anisotropy class yields 
Nc — 9.44121. The equivalent of (|6j becomes dig — 
eg + 0.549(A^ - Nc)g^ + A7.6g^, leading to d^^ w 4 - 
0.00158(A^ - Nc)^. Although it yields di^N = 3) w 3.93 
and no QLRO phase in d = 3, naive extrapolation should be 
taken with caution given the high value of Nc- Numerical val- 
ues for die are changed for a 7^ ^ 0, but the scenario is robust 
for7Q < 7c 133]. 

We now discuss the FRG flow-equations for A^ large. From 
a truncated exact RG Tarjus and Tissier (TT) |26] found: that 
the linear cusp of the F/D fixed point for d > 4 vanishes for 
N > N*{d), i.e. i?"'(0+) = 0; and that the non-analyticity 
becomes weaker as A^ increases (as |0|" with n ^ N). An- 
alytical study of the derivatives of (|2ji confirms the existence 
of this peculiar FP to two loop and predicts A^* (d, 2p) beyond 
which the set of {i?^^'^' (0)} for A: < p admits a stable FP, with 
^(2fe-i) (0+) = for fc < p and i?(2fc-i) (0+) ^ for k > p. 
We find: 

iV*(d) = Ar*(d,4) = 18 + 49e/5+... (9) 

which yields a slope roughly twice the one of Fig. 1 of J2^. 
This remarkable FP raises some puzzles. Although weaker 
than a cusp its non-analyticity should imply some (weaker) 
metastability in the system. It is thus unclear whether DR is 
fully restored: to prove it one should rule out feedback from 
anomalous higher-loop terms in exponents or the /3-function. 
Finally, one also wonders about its basin of attraction. As 
shown in Fig. 13 flie FRG flow for R""{0) is still to large val- 
ues if its bare value is large enough, indicating some tendency 
to glassy behaviour 

To explore these effects we now study the F/D phase tran- 
sition at large A^ and d > 4. We obtain, both at large A^ and 
fixed d (extending Ref. 1E3I1 ). and to one loop, the flow equa- 
tion for the rescaled R{z = cos(/)) = NR{(t))/\e\: 

diR= -R + 2R[R- R[R'z + ^R'^ ^0 . (10) 

We denote y{z) = R'{z), yo = ^'(1) = -NR"{0)l\e\ and 
r4 = NR""{0)/\e\. There are two analytic FPs R{z) = 
z — 1/2 and R(z) — z'^/2, corresponding both to j/o = 1 
and to r4 = 1 and ~ A respectively. This agrees with the 
flow of thed erivatives for analytic R{(j))'- diyo = 2/0(2/0 ~ 1). 
and at yo = 1: dir4 = ■i(r4 — l)(r4 — 4). The first FP is the 
large-A^ limit of the TT fixed point, the second is repulsive 
and divides the region where r4 — > cx3 (non-analytic i?((/')) in 
a finite RG time Ic (Larkin scale). For 2/0 > 1, we find a fam- 
ily of NA fixed points with a linear cusp, parameterized by an 
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R""{0) 




FIG. 3: Flow for i?""(0)/|e| for d = 4 - e > 4 as afunction of iV. 
The two branches behave as and 4/7V at large TV. 



integer n>2, s.t. yo ^ n/{n~ 1), z = y - {ya - l){y/ yo)"- 
The solutions with n (i.e. z{y)) odd correspond to random 
anisotropy iR'{(f)) = R'{(f) + tt)). The n = 2 RF fixed point 
is = 2cos((/)) + ^ sin^((/)/2) - |. To elucidate their 
role, we obtained the exact solution for the flow both below 
Ic, i.e. z = ^ + (j/o — (*^(^) parameterizes the bare 
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disorder, $(1) = 0), and above Ic, with an anomalous flow for 
yo. Matching at Ic yields the critical manifold for RF disorder, 
defined from the conditions that ^'{w) = ^{w)/w = 1 has a 
root < u) < 1. It is different from the naive DR condition 
yo — 1, valid for small r4. The n = 2 FP corresponds to bare 
disorder such that the root w = 0. Hence it is multicritical 
js^l- Generic initial conditions within the critical F/P mani- 
fold flow back to the TT FP i.e. the linear cusp decreases to 
zero ii^ . This however occurs only at an infinite scale, hence 
we expect a long crossover within a glassy region, character- 
ized by a cusp, and metastability on finite scales f3^. The 
large-A^ limit here is subtle. Taking iV ^ cx) at fixed volume 
on a bare model with Ro{z) — z yields only the analytic FP, 
equivalent to a replica-symmetric saddle point. Higher mono- 
mials z'P are generated in perturbation theory, at higher order 
in 1 /N. Thus, for N large but fixed and infinite size, one must 
first coarse grain to generate a non-trivial function Ro{z), be- 
fore taking the limit of iV ^ oo. 

In conclusion we obtained the 2-loop FRG functions for 
the random field and anisotropy cr-models. We found a new 
fixed point and a scenario for the decrease of the lower critical 
dimension. This rules out the scenario left open at one loop 
that the bifurcation close to d = 4 simply occurs within the 
(quasi-) ordered phase. 
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